KA-TP-09-2009 



Next-to-leading order QCD corrections to 
W + W + jj and W~W~jj production via weak-boson fusion 

B. Jager 

Institut fur Theoretische Physik und Astrophysik, 
Universitdt Wiirzburg, 
97074 Wiirzburg, Germany 

C. Oleari 

Universitd di Milano-Bicocca and INFN Sezione di Milano-Bicocca, 

20126 Milano, Italy 

D. Zeppenfeld 

Institut fur Theoretische Physik, 
Universitdt Karlsruhe, KIT, 
76128 Karlsruhe, Germany 

Abstract 

We present a next-to-leading order QCD calculation for e + v e fi + v tl jj and e~v e n~v IJi jj production 
via weak boson fusion at a hadron collider in the form of a fully-flexible parton-level Monte Carlo 
program, which allows for the calculation of experimentally accessible observables within realistic 
selection cuts. The QCD corrections to the integrated cross sections are found to be modest, while 
the shapes of some kinematical distributions change appreciably compared to leading order. The 
residual scale uncertainties of the next-to-leading order results are at the few-percent level. 
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I. INTRODUCTION 



Weak-boson fusion (WBF) reactions have been identified as a promising means for gaining 
a thorough understanding of the electroweak gauge boson sector of the Standard Model 
(SM) and extensions thereof. WBF is considered as a possible discovery channel for the 
Higgs boson and a suitable tool for a later determination of its couplings to 



gauge bosons and fermions 



0, 0,3 



and of its CP properties 



0,Q. 



As WBF reactions 



are sensitive to the mechanism of electroweak symmetry breaking per se, they could help 



to spot si 
sector [11 



natures of physics 



12. 



13 



14 



15 



16 



jeyond the SM such as strong interactions in the electroweak 



17 



181 ]. Furthermore, WBF processes with like-sign dilepton 



final states are an important background to various new physics scenarios. 

In order to unambiguously identify new physics effects, precise predictions for signal and 
background processes are crucial, including estimates of the theoretical uncertainties. In the 
context of perturbation theory, the required accuracy and a first assessment of the uncertain- 
ties associated with the truncation of the perturbative expansion can only be achieved by the 
calculation of next-to-leading order (NLO) QCD corrections for experimentally accessible 
observables. NLO-QCD corrections have therefore been provided for several WBF-induced 



processes within the SM, including electroweak H 



and W ± Zjj production at hadron colliders 



19 



20 



7.7 



21 



HirLW 



22 



23 



24 



rL Z33, W+W-jj, ZZjj, 



251 ] . QCD corrections have 



also been determined for selected WBF processes in new physics scenarios, such as [26l . |27| . 
More recently, the electroweak corrections [28J] to the Higgs boson signal in WBF have been 
computed, together with parts of the two-loop QCD corrections 29j and the interference 



with QCD-induced Hjj production 



30, 



31] . Supersymmetric corrections to pp — > Hjj have 



been presented in the Minimal Supersymmetric Standard Model (MSSM) 32]. 

In this work, we focus on e + v e ^i + v^jj production via WBF at a hadron collider and 
present results for the CERN Large Hadron Collider (LHC). We develop a flexible Monte 
Carlo program that allows for the calculation of cross sections and kinematical distributions 
within arbitrary selection cuts at NLO-QCD accuracy. We discuss the most important theo- 
retical uncertainties and quantify the impact of radiative corrections on a few representative 
distributions. 

Results for the e _ ^ e /i~i^jj channel can simply be obtained by charge conjugation and 
parity reversal. For example, if one is interested in pp —>■ e'v^v^jj at the LHC, it is 
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sufficient to run our code for the process pp — > e + v e jj, + v^jj with a centre-of-mass energy of 
14 TeV, treat the final-state charged leptons as if they were negatively charged, and reverse 
momentum directions when considering parity-odd distributions. For this reason, in the rest 
of the paper we will discuss only the e + u e fi + i , fl jj process. 

We start with a brief overview of the technical pre-requisites of the calculation in Sec. HT1 
Section [TTI1 contains our numerical results. We conclude in Sec. [TV] 



II. TECHNICAL PRE-REQUISITES 

The calculation of NLO-QCD corrections to pp — ► e + h , e fi + u ll jj via WBF proceeds along 
the same lines as our earlier work on WBF W + W~ j.j, ZZjj, and W^Zjj production in pp 



collisions. The methods developed in Refs. 



23 



24 



251 ] can therefore be straightforwardly 



adapted and only need a brief recollection here. 

At order a 6 , WBF e + iy e ^i + u fl production in association with two jets proceeds via the 
scattering of two (anti-) quarks by t-channel exchange of a weak gauge boson with subsequent 
emission of two W + bosons, which in turn decay leptonically. Non-resonant diagrams where 
leptons are produced via weak interactions in the t-channel also have been included. For 
brevity, we will refer to pp — > e + u e fi + u^jj as "WBF W + W + jj" production in the following, 
even though the electroweak production process includes non-resonant diagrams that do not 
stem from a W — > iv decay. 

For each partonic sub-process, the 93 contributing Feynman diagrams can be grouped 
into four topologies, which are depicted for the representative uc — > dse + v e fi + i> ll mode in 
Fig. [TJ The first two topologies correspond to the emission of two weak bosons from the 
same (a) or different (b) quark lines with the W + decaying leptonically. Topologies (c) 
and (d) comprise the so-called "leptonic tensors" Lww and Tyw,e, which describe the tree- 
level amplitudes of the sub-processes W + W + — > e + z/ e /i + z/^ and VW + — > e + v e , respectively. 
In each case, V stands for a virtual photon or Z boson. Not shown are diagrams analogous 
to (a) and (d) with W + emission from the lower quark line rather than the upper one and 
graphs which are obtained by permuting the final-state leptons. 

We disregard graphs with a weak-boson exchange in the s-channel with subsequent de- 
cay into a pair of jets, as they are strongly suppressed in the regions where WBF can 
be observed experimentally. For sub-processes with identical-flavor combinations such as 
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FIG. 1: Feynman-graph topologies contributing to the Born process uc — > <ise + f e /x + f M . Not shown 
are diagrams analogous to (a) and (d) with W + emission from the lower quark line. V denotes a 
Z boson or a photon. 

uu — > dde + h> e fi + v t j i , in addition to the t-channel contributions discussed above also w-channel 
diagrams and their interference with the ^-channel graphs arise. While we do take into ac- 
count both t- and ^-channel contributions, their interference cross section is kinematically 
strongly suppressed {28] and therefore not considered in the following. 

The relevant diagrams are combined in an efficient way, avoiding multiple evaluation 
of recurring building blocks, and numerically evaluated using the amplitude techniques of 



Refs. 



33 



34|. 



For the NLO-QCD calculation, real emission and virtual corrections have to be consid- 
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ered. Infrared singularities are handled via the subtraction formalism in the form proposed 



by Catani and Seymour 



35]. 



The real-emission contributions are obtained by attaching one extra gluon to the tree- 
level diagrams sketched above. In addition to (anti-) quark initiated sub-processes such as 
uc — > dsge + v e fi + i> fl , contributions with a gluon in the initial state (e.g. ug — > <isce + z/ e /i + z/ M ) 
emerge. 

The virtual corrections comprise the interference of one-loop diagrams with the Born am- 
plitude. They are calculated in d = 4 — 2e dimensions in the dimensional reduction scheme. 
Due to color conservation, contributions from graphs with the virtual gluon being exchanged 
between the upper and the lower quark line vanish at order a s , within our approximations. 
Only self-energy, triangle, box and pentagon corrections to either the upper or the lower 
quark line have to be considered. The singularities stemming from infrared divergent con- 
figurations are calculated analytically and canceled by the respective poles in the integrated 
counter-terms of the dipole subtraction approach. The finite pieces of the loop diagrams 
can be calculated in d — 4 dimensions. The emerging two-, three-, and four-point tensor 
integrals are evaluated numerically by a Passarino-Veltman reduction procedure. In order to 



avoid numerical instabilities, for the computation of the pentagon integra 
reduction scheme of Refs. 



36 



25 



s we resort to the 
. We monitor 



371 ] , which has already been employed in 
the numerical stability by checking Ward identities at every phase-space point and find that 
the fraction of events which violate electroweak gauge invariance by more than 10% is at the 
permille level. The respective phase-space points are disregarded for the calculation of the 
finite parts of the pentagon contributions and the remaining pentagon part is compensated 
for this loss. As the pentagon contributions amount to only about 3%o of the full NLO result, 
the error induced by this approximation on cross sections and distributions is negligible. 

In order to ensure the reliability of our calculation, in addition to the stability tests for 
the pentagon contributions, we performed the following checks: 

• We compared our tree-level and real-emission amplitudes for WBF e + h , e fi + i , fl jj scat- 
tering and crossing-related processes to the corresponding expressions provided by the 



automatized matrix element generator MadGraph 39] and found full agreement within 
the numerical accuracy of our program. 

We compared our integrated cross section for pp — > e + v e ^, + v^jj via WBF within the 
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approximations discussed above to the full cross section provided by MadEvent [39] 
within typical WBF cuts. Our predictions fully agree with the corresponding MadE- 
vent results within the statistical errors of the two calculations. Agreement between 
the two programs is also found for representative distributions. 



III. NUMERICAL RESULTS 



The cross-section contributions discussed above have been implemented in a fully flexible 
Monte-Carlo program, structured analogous to the VBFNLO code [40]. In order to facilitate a 
comparison of the general features of this production channel, similar settings and selection 



cuts as in Refs. 



23 



24 



251 ] are employed. We use the CTEQ6M parton distributions 41] 



with a s (mz) = 0.118 at NLO, and the CTEQ6L1 set at LO. Since in our calculation quark 
masses are neglected, we entirely disregard contributions from external b and t quarks. As 
electroweak input parameters, we have chosen irtz = 91.188 GeV, mw = 80.423 GeV, and 
Gp = 1.166 x 10~ 5 / GeV 2 . The other parameters, oqed and sin 2 ^, are computed thereof 



via tree-level electroweak re 
use the algorithm 



42 



43. 



ations. For the reconstruction of jets from final-state partons we 
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451 ] with resolution parameter D = 0.7. All our calculations 



are performed for a pp collider with center-of-mass energy of \fs = 14 TeV. 

In order to clearly separate the WBF signal from various QCD backgrounds, we employ 
a set of dedicated selection cuts. We require at least two hard jets with 



p Tj >20GeV, |y i |<4.5. 



(1) 



where pTj denotes the transverse momentum and yj the rapidity of a jet j with the latter 
being reconstructed from massless partons of pseudo-rapidity \rjj\ < 5. The two jets of 
highest transverse momentum are referred to as "tagging jets" . We impose a large rapidity 
separation between the two tagging jets, 

A Vjj = IVji-vM >4, (2) 

and furthermore demand that they be located in opposite hemispheres of the detector, 

Vji x Vj2 < , (3) 

with an invariant mass 

Mjj > 600 GeV . (4) 



For our phenomenological analysis we focus on the e + u e fi + u fl leptonic final state. Multi- 
plying our predictions by a factor of two, results for the four-lepton final state with any com- 
bination of electrons and/or muons and the associated neutrinos (i.e., e + z/ e /i + z/ At , e + v e e + v ei 
and fi + u fJt fi + u t j) can be obtained, apart from numerically small identical-lepton interference 
effects. 

To ensure that the charged leptons are well observable, we impose the lepton cuts 

Pr/>20GeV, M<2.5, (5) 
AR je >0A, AR U >01, (6) 

where ARj£ and ARu denote the jet-lepton and lepton-lepton separation in the rapidity- 
azimuthal angle plane. In addition, the charged leptons are required to fall between the two 
tagging jets in rapidity, 

yj,min ^ T}£ Vj^max- (7) 

As discussed in some detail in Ref. 



251 ]. a proper choice of factorization and renormal- 
ization scales, \ip and [ir, can help to minimize the impact of higher order corrections on 
cross sections and distributions in WBF reactions. In the W^Zjj case, the momentum 
transfer Q between an incoming and an outgoing parton along a fermion line has been found 
to be more suitable than a constant mass scale. A similar behavior can be observed in the 
W + W+jj CclSG, clS illustrated by Fig. [2J where we show the total cross section within the 
cuts of Eqs. ([I])-©, <r cuts , as function of fip and //r, which are taken as multiples of the 
scale parameter /x , 

= £,f^o , V-r = £,r^o ■ (8) 

We distinguish between /x = m w and fio = Q and vary £p and in the range 0.1 to 10. At 
LO, the variation of the cross section with £p resembles the scale dependence of the quark 
distribution functions, f q (x,fip), at relatively large values of x. Beyond LO, cr cuts barely 
depends on //p. The renormalization scale enters only at NLO via the strong coupling 
a s (fj,R). Due to the small size of the radiative corrections, the resulting /i# dependence of 
the NLO cross sections is marginal. In Tab. [I we list the results for er cuts at LO and NLO 
for the two scale choices (jlq = mw and [Iq = Q with £r = £f = 1 and the settings defined 
above. For reference, we also quote the respective results obtained with the new MSTW08 
parton distribution functions of Ref. j^. It is interesting to note that these results differ 
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FIG. 2: Dependence of the total WBF pp — > e + z/ e /i + z/ At jj cross section at the LHC on the factoriza- 
tion and renormalization scales for the two different choices of ^o : Mo = m w and Mo = Q- The NLO 
curves show <r cuts as a function of the scale parameter £ for three different cases: fiR = ^f = £mo 
(solid red), = ?Mo and fiR = (dot-dashed blue), hr = ^fio and \if = Mo (dashed green). The 
LO cross sections depend only on \ip (dotted black). 



PDF set 


°lo s (mo = m w ) 


ct nlo(Mo = mw) 


alt s (V0 = Q) 


< s o(Mo = Q) 


CTEQ6 


1.478 fb 


1.404 fb 


1.355 fb 


1.402 fb 


MSTW08 


1.409 fb 


1.372 fb 


1.287 fb 


1.369 fb 



TABLE I: Cross sections for pp — > e + v e n + v^jj via WBF at LO and NLO within the cuts of 
Eqs. CEE]) for the scale choices fip = Hr = mw and fJ>F = fJ-R = Q and two different sets of parton 
distribution functions. The statistical errors of the quoted results are at the sub-permille level and 
therefore not given explicitly. 

from the corresponding cross sections obtained with the CTEQ6 sets by about 5% at LO 
and 2 to 3% at NLO. In the following we will stick to mf = M-R = Q and the CTEQ6 parton 
distributions. 

A very characteristic feature of WBF-induced processes is the distribution of the tagging 
jets, which are located in the far-forward and -backward regions of the detector. Due to the 
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FIG. 3: Rapidity distribution of the jet with the highest (a) and with the lowest (b) value of the 
absolute value of the rapidity \y^\ in WBF pp — > e + v e ^ + v^jj production at the LHC at LO (dashed 
black line) and NLO (solid red line). 

color-singlet nature of the t-channel weak boson exchange, the central-rapidity range of the 
detector exhibits little jet activity. As this feature is an important tool for the suppression 
of QCD backgrounds which typically exhibit rich hadronic activity at low rapidities [4], it is 
important to estimate the impact of NLO-QCD corrections on the rapidity distributions of 
the observed jets. To this end, in Fig. [3] (a) we display the distribution of the jet with the 
highest absolute value of rapidity at LO and NLO. The respective results for the jet with 
the lowest absolute value of rapidity are shown in Fig. [3] (b). While at LO each final-state 
parton corresponds to a tagging jet, at NLO two or three jets may arise due to the extra 
parton in the real-emission contributions. Because of this extra parton, the probability to 
find a jet of very high or very low rapidity is larger at NLO than at LO. 

In Fig. HI the invariant mass distribution of the two tagging jets is shown for two choices 
of /io together with the dynamical K factor, which is defined via 

_ da NLO /dM j:j 



K{M. 



33 J 



(9) 



^o"lo / dMjj 

The Mjj distribution was scrutinized in Ref. [25J to demonstrate the impact of the choice 
of factorization and renormalization scales on the perturbative stability of calculations for 
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FIG. 4: Invariant mass distribution of the tagging jets in WBF pp — ► e + v e fi + u fJi jj production at 
the LHC for two different choices of /^o [panels (a) and (b)] at LO (dashed) and NLO (solid). Their 
ratios, the K factors as defined in Eq. Q, are displayed for /iq = mw in panel (c) and for Ho = Q 
in panel (d). 

WBF processes. Using the example of WBF W + Zjj production it was shown that small 
NLO-QCD corrections to da/dMjj are obtained, if /j,p and are identified with Q. A 
similar behavior is observed for the W + W + jj case, with a K factor close to one over the 
entire range of Mjj for /i = Q, but pronounced shape distortions at LO for /i = m w . 
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IV. SUMMARY AND CONCLUSIONS 



In this work, we have presented a NLO-QCD calculation for WBF W + W + jj and 
W~W~jj production at the LHC, which takes leptonic decays of the weak gauge bosons 
fully into account. We have developed a flexible Monte-Carlo program that allows for the 
calculation of cross sections and distributions within typical WBF cuts. The QCD correc- 
tions to the integrated cross section are modest, changing the leading-order results by less 
than about 7%. The residual scale uncertainties of the NLO predictions are at the 2.5% 
level, indicating that the perturbative calculation is under excellent control. In accordance 
with our earlier work on WBF W ± Zjj production we found that size and shape of the NLO 
results can be best approximated by the corresponding LO expressions, if the momentum 
transfer, Q, of the t-channel electroweak boson is chosen as factorization scale. 
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